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ABSTRACT 


The failure point of a unidirectional composite subjected to longi- 
tudinal shear loading is calculated. The method of analysis is based on 
the finite element technique, the incremental plasticity relations of 
Prandtl-Reuss, and the von Mises yield criterion. The failure point for 
a boron-aluminum composite and a boron-epoxy composite were eorouced to 


determine the effect of the matrix material on the composite. 
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I INTRODUCTION 


The need for material possessing a higher strength to weight ratio 
has become abundantly clear to modern industry. Existing materials 
satisfy the present high strength-low weight requirements (high per- 
formance), but composite material technology, still in its infancy 
promises to raise the standards of high performance even higher. Com- 
posite material technology is particularly interesting to the aerospace 
industry which is constantly striving to improve the performance of aero- 
nautical and aerospace vehicles. Therefore a possible solution to the 
aerospace industries’ quest for higher performance materials is afforded 
by composite materials. A composite material combines the high strength 
and high modulus of one material and the low-strength and low modulus 
of another material (or materials) to yield a material exhibiting the 
necessary high strength-low weight ratio Hes ined by the aerospace 
industry. 

A particularly useful composite is the two material fiber reinforced 
composite which consists of a filament embedded in a compatible matrix 
material. The filament possesses the high strength and high modulus 
properties required while the matrix, having a low strength and low mod- 
ulus adds the ductility necessary in a high strength low weight ratio 
composite. The matrix also provides the form of the composite required 
for a particular structural use. A composite plate structure usually 
consists of many layers of composite lamina with the adjacent lamina 


having different filament orientations. 


The composite structure is analyzed considering each lamina as an 
anisotropic homogeneous material, which results in the calculation of 
macrostresses. If a lamina is examined on an individual point by point 
basis, and considering the difference in the mechanical properties of 
the filament and the matrix, the stresses calculated are ihemseanle. 
The microscopic stresses are important in determining at which point in 
the composite that yield or failure occur. 

To determine the yield or failure point of the composite, tensile 
tests are performed on the matrix and filament materials until the 
elastic and plastic limits are reached. The value of the limits, 
namely the yield stress and failure stress are microscopic stresses. 

To accurately determine the behavior of a composite material the macro- 
scopic stress state must be related to the microscopic stress state. 
This may be accomplished using various analytical approaches. In this 
Study the finite element technique is employed. 

Several authors have considered the macrostress analysis of uni- 
directionally reinforced composites. Bloom and Wilson [4] have calculated 
the elastic microstresses for longitudinal loading. Longitudinal shear 
loading analysis was obtained by Adams and Doner [1], Baker and Foye [3] 
and Tsai, Adams, and Doner [15]. The transverse normal loading problem 
was considered by Foye [8], as well as Adams and Doner [2]. Lin, 
Salinas, and Ito [11] employed the finite element technique to solve the 
problem of combined loading, calculating an initial yield surface and 
later extending their work into the plastic region [10]. 

The present investigation considers the effect of longitudinal 
Shear loading on a boron-aluminum composite and a boron-epoxy composite. 


The objective of this analysis is to determine the failure loads of the 
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two composites and to compare the effect that the two different matrix 


materials have on the composite's properties. 


im 





Mie VesCRIRMONBOr PLASTICITY 


This investigation considers the elastic-plastic analysis of a uni- 
directional reinforced composite material subjected to longitudinal 
shear loading. Elastic theory being a basic subject and plastic theory 
being a more difficult, less studied topic will be considered here and 
amplified in section three, while the elastic theory will only be 
briefly discussed in section three. This section describes some of the 
assumptions of plasticity while also describing the yield criterion 
employed in this analysis, the Prandtl-Reuss equations and the universal 


stress-strain law assumption. 


“A. YIELD CRITERION 

Plasticity may be defined as that property that allows a materia! 
to deform continuously and permanently without rupture during the 
application of stresses that exceed those producing yield. The final 
distortion does not only depend on the final state of stress, but upon 
the sequence of stress states producing that final state. This is in 
contrast to elastic deformation which depends only on the final stress 
state. 

An important assumption of the theory of plasticity of metals, 
based on experimental observations, is that the volume of material 
remains constant under plastic deformation, that is, a hydrostatic 
pressure does not cause yielding. This implies plastic strain 1S aSsoci- 
ated with distortional energy only and hence the hydrostatic component 
of a state of stress does not influence the point at which yielding 


occurs. What 1s needed then, is some criterion, a so called yield 
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criterion, that will predict yielding under a state of stress, given 
only the yield stress determined from a simple uniaxial tension test. 

For a homogeneous, isotropic material (same properties at all points 
and all directions), a criterion should not allow hydrostatic stress to 
influence yielding. Yield must also be independent of the directions 
of the axes chosen to define the system. Therefore since the material 
is isotropic, yielding will be related only to the intensity of stress 
and must be a function of the invariant of the stress tensor. The 
invariants are the coefficients of the cubic equation, which is associ- 
ated with determining the principal stresses of the stress state. If 
yielding is unaffected by hydrostatic stress (P) then it depends only 
on the deviatoric stress. This means that the criterion must be a func- 
tion of the deviatoric invariants. 

Before proceeding with the yield criterion discussion, some terms 
need explaining. Hydrostatic stress, P, is a stress that acts equally 
in all directions and will cause only a change of volume, recoverable 
on removal of the stress. This is related to the dilatation, A which 
is defined as the change of volume per unit volume. eager tone 
stress components, Cage are obtained as the difference of the dilatation 


from the total stress tensor, that is 


Ae eS ope 
O5 644? O55 
or 
a) aaa aid a 855 P (2-1) 
ae . 
where P = Te and Fe 1s the Kronecker delta. 


Upon expansion, equation (2-1) gives 


ae Cera (2-la) 


» 5 a 2-1 
Syy = Syy"P (allo) 


ils 





oe = Gog (2-Ic) 
ey = Sivay (2-1d) 
Ov5 = OV, (2-le) 
O5y Coy (2-1f ) 
where the hydrostatic stress, P, may be expanded to 
°F 3 
ha (o, to, to, 71/3 (2-2) 


where the Einstein summation convention is followed, where a repeated 
subscript means summation with respect to the subscript. In this study 
index notation as well as cartesian notation will be exchanged freely 
instead of the indices having the usual numerical values. The devi- 


atoric stress invariants are also expressed in this manner. 


| “4. —_ cet 
J “= ore TF yy Sy" =0 (2 3a) 
4x 2g. 7! e222? , - 22 15 _2h\ 
us 2744715 BLY TOV 7 aay iy OTe )] (2-3b) 
J3 "50% =[o. 7 *+20.70. 76.7 


1 ko ki XX on TAFE SESE TON 


=o eye =e FO ae “975 ie _ (2-3) 
Various yield criterion exist that were verified experimentally and 
that are applicable to the problen. The von Mises yield criterion is 
the simplest and most convenient to use. It may be expressed as a 
function of the deviatoric stress invariants 
F(J,~)=0 
or specifically we will use 
Jo* = k? (2-4) 
where k is a constant. When Jo“<k* there is elastic behavior, that is 


there is no change in the plastic strain. When Jo“=k? there can be 


elastic or plastic strain depending upon the rate of change of Jo”. 
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Other derivable forms of J»~ follow in terms of the principal stresses 


Ols O95 O3- 


205 °=01 “240972403 72 2k2 (2-5a ) 
or 
239°=4 (1-62 )?# (69-03) 2# (03-01 )21=2k? (2-5b) 
also 
Aull = 2 _ 2 s 2 2 2 2a 912 2 
205 Sal Coos ov) ING coe) A Gia oye) ]+2Loy toy, oe J=2k2 (2-6) 


The constant k may easily be solved for in terms of Y, the yield 


stress in simple tension, by substituting into equation (2-5b), o,=Y 


and O5=03=0 
01 *+o,7=6k2 
or 
2Y2=6k2 and Kae (2-7) 
V3. 
Now consider the meaning of k in pure shear. Pure shear is equiv- 
alent to the stress state roe ie o,>0 and o3=-0, with oo=9. Substituting 


into equation (2-5b) again 
2 

Cae nay 

and (2-8) 


Mone 
Tmax V3 


Hence k is the shear stress, at which a point in pure shear, will yield. 


Ca iar lp) 
) 2+ ( Ne) 6k 


The von Mises criterion implies therefore, that yielding is not 
dependent upon any particular stress component, equation (2-5b) shows 
that it depends upon a function of distortion since the differences of 
principal stresses are proportional to the maximum shear stresses. 

It is convenient to state yet another variation of the yield cri- 
terion. From equations (2-5b) and (2-7), we obtain 


aL 
Y¥=={(0q-09)2+(09-03)2#(03-0,)2} = Go (2-9) 





The scaler quantity on the right hand side may be defined as an 
equivalent stress, a; so that when c=Y, yielding occurs. Thus, any 
value of o less than Y will give rise to elastic states of stress only. 
When s=Y there may be either elastic and or plastic components, depending 


upon the sign of do. 


B. FLOW RULE 

Having selected a yield criterion it 1s necessary to establish a 
relationship between the stress components and the corresponding deform- 
ations. This is usually referred to as the flow rule. 

St. Venant proposed that the directions of the increments of prin- 
cipal plastic strains correspond to the directions of the principal 
stresses. Later Levy proposed a relationship to connect the plastic 


strain increment to the stress in the form 





X 
d)s=-——— FOO Oe FE (2-10) 


where the superscript P denotes plastic. ; 
The factor dr establishes the relation between quantities of the 
same differential order. From the first and third equations of (2-10), 


P Po 
de de =) ( 


; P nee 7 7" 
yy a ae de -de,, dro, aes). (2-11) 


Adding the two equations of (2-11) and recalling that the volume change 


a P_,_P ee 
due to plastic strain is zero, that is, dV sdey tds ide | 0, or, 
P a Pp | 
de tds, ec. we obtain 
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Ws 
eedy Gg 4 


d : 
XK EXX 2 


“yy *8z7) 

6 Equations (2-12) 
De 

de Oro, 


Similar equations may be obtained for the other components by cyclic 


permutation of (2-12). 


These equations may also be written in deviatoric form, using 
2 
0 ( ) =2f 


~ 1 
XX OXX (oyyto JJ 


Ge =O. to COC = 
XX YY FG h NOG 
: : 3 
Ome 
XY XY 


and comparing with equation (2-12), we obtain 


Be Y 
ae x 
lg laa : 
de vy dhe 
or, in general form 
a , pie 
de; ; dha 4 (2-13) 


In the above equations, a rigid plastic material is assumed because 
no mention is made of any elastic strains. For a complete statement, 


the equation should be written : 


; P = a _ 
de 54 dross (2-14) 
where the superscript P denotes the increment of plastic strain. The 


total strain increment becomes 


. P 
des des, + de, . 


(2-15) 
where ort is the elastic strain increment. Now, according to the 
present argument, (2-14) refers to increments of the strains as we 
follow the deformation path of the material; we must therefore add the 


corresponding elastic strain increments - that is, although the plastic 


Strain increments are, by definition, proportional to the current values 


17 








of the deviatoric stress o.7, the contribution of the elastic strain will 


1J’ 
be those corresponding to the change in the stress state as we go from 
two closely neighboring states of stress. 

The complete statement of the stress-strain relation becomes 


de. area 2-16) 
Equation (2-16) is the Prandtl-Reuss equation for the deviatoric strain 
increment. The equation derives its name from Prandtl who developed the 
elastic and plastic strain it and a) equations for plane stress and 
Reuss who generalized the equation, but it was Lévy who originally de- 
veloped the plastic strain equation. 


The rule of plastic flow in terms of principal stresses and strains 


also gives 


tes * = 3 = qy (2a) 











C. THE UNIVERSAL STRESS-STRAIN LAW 

The Universal stress-strain law is an assumption that the plastic 
uniaxial stress-strain relation Coe f(de") is the sane relation be- 
tween equivalent stress and strain, that is 6 = £(de?) for the general 
stress state. proto ing, cie universal stress-strain law equation (2-17) 
may be expressed in terms of equivalent (effective) strain and stress as 


-c¢& 
cas (2-18) 


The constant C depends upon the definitions of equivalent strain and 
stress. The equivalent strain increment may be expressed 
al | f2 2 2 2 
de =|? {(de,?) + (dep) a (de!) 2 
i 


P ix 
C43 de; 2 | DET) 


4 


ie 
3 
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(2 vo 
The numerical value of 3 has been chosen so that in uniaxial stress, 


See 


Pp Pas P Pa P 
(where dev = eee ; de.) we have de dey. The assumption of a 
universal stress-strain law then gives a value of C=°/5. The actual 


proof of the numerical value of C may be found in Appendix B. 


ie 





TTI DERIVATION OF THE STRESS-STRAIN MATRIX 


The purpose of this section is to give an explicit expression for 
ro’ the plastic stress-strain matrix for the von Mises material. The 
expression obtained takes .a form that can be accommodated to the finite 
element analysis. | 

The following derivations follow closely those presented by Yamada, 
Yoshimura and Sakurai [16]. The well known equations of elasticity, 
commonly called Hooke's law may be expressed in matrix form for an 


isotropic material as 
































a l-v v v a 

o" meer ee) Wes a 
s v l-v v z 

YY Tey 1 yy 
He F v v ]-v : 
oa jay | U2 Te2v 1-2 aes 
Oy 0 0 0 CExy 
C7 0 0 0 cev, 
o 0 0 0 ce 

Zx le \e zx | 


(3-1) 
where E is Young's modulus of elasticity and v is Poisson's ratio. The 
matrix equation may be expressed as 

(o}=E[D =] fe} or since E=2(1+v)/G 
where G is the shear modulus of elasticity 
{o}=2(14+v) (0 ]{c} (3-2) 
where roe] is referred to as the elastic stress-strain matrix, {co} as 


the column matrix of stress and {e} as the column matrix of strain. 
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Hooke's law for the longitudinal shear loading case (i.e., all strains 


equal to zero except Eve and Hr) greatly simpifies equation (3-1) to 


| a | 
Oo. = § Ce 
y2 li 2 (3-3) 
0 


— 


Z| 172 
It is interesting to note that in the elastic range xz and yz behavior 
are independent. This is illustrated by equation (3-3) where the 
equations for Os and oS uncouple. 
The inverse of Hooke's law, equation (3-1), may be expressed in index 
notation as 
‘eure o 
eye es ao (qeey) abs (3-4) 
ers "J 3E 
where E45 are the components of the strain tensor. The O45 are the 


components of deviatoric stress which are defined as 


ee 


= _) 
i pe ae: (3-5) 
and explicitly expressed by equation (2-1). 
Before continuing with the derivation of the stress-strain matrix 


and deviatoric strain components, E44 must be explicitly stated, 


~ 


Sccimimea oial (3-6a ) 
Fyy = Fyy® sic 
2 ie Cea (3-6c) 
Exy = Exy (3-6d ) 
Ev? = Evy (3-6e) 
Bay ~ Ezx cao 
where the dilatation strain is 
g = 62 35, EN (3-69 ) 


Equations (3-6) will help better understand the Prandtl-Reuss equations 


for the deviatoric strain increment expressed in index notation as 


2] 








— r alos 5 
2G 
where according to equation (2-18) 
‘dy = c-2E 
6 
Letting H- ale, dx becomes 
de? dé 
dx = Cl (3-8) 
oH~ 
and allowing the constant C to have the value three halves 
eR (3-9) 
oH~ 


The equivalent plastic strain increment dle and the equivalent stress 


c are equal to 


oe 1 
Pp Pp Pi? 
de = (4 dej, de. Ne (3-10) 
1 
APS. fee - 
do 
H- = —5 shown in equation (3-8) is actually the slope of the stress 
de 


versus plastic strain curve in the non-linear plastic region (increment 
of equivalent stress (da)/increment of equivalent plastic strain me 
This is shown by figure 1, where for the plastic problem the transition 
must be made from the total equivalent stress-strain curve on this 
universal stress-strain curve to the equivalent stress-strain curve for 
the non-linear range. 

For the derivation of the stress-strain matrix for this elastic- 
plastic analysis of the longitudinal shear problem the equivalent stress 
becomes 


2 2 
2 = oe e 
S {3(o,5 ws, (3-12) 


ee 


and the plastic strain increment dee becomes 


as 1 
ae, 2 2 2 
de = {3 (de ps tde Po )} (3-13) 
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FIGURE 1. SLOPE OF THE NON-LINEAR STRESS-STRAIN CURVE. 


ooo 


The von Mises yield criterion is employed in the plastic analysis 


of the problem. It is expressed here in terms of the deviatoric stress 


invariant 


205° B™BO.: 0. - oan (3-14) 


ij 45 


wilh 


and in differential form 


“do.t = <add 35 
eae 43 a 30d0 ( = ) 
Solving equation (3-9) for do and substituting into equation (3-15) we 
have 

| a c Sao = = 
Algebraically eliminating do. 3 from equations (3-7) and (3-16) gives 


, 2 , * ae) ’ = 
eGo. : (do, 5 05 3dX) g0“H"d) (3-17) 


Using equation (3-11) and solving for dd 


. (3-18) 
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For convenience let 


S = £52(1+H~/36) 


thus 
dri = O45 des; 
(3-19) 
S 
The term O45 Ges: may be simplified to O45 Cea The proof is not very 


long and will be presented here rather than detailing it to an appendix. 


045 qe} 5 = oe dey + Gy dev, * cme des 


#2(0 xy dey + °yz dev, * Lee, 
expanding 


| Oo ar 2e css 
= 72 ee XX 


Oo. +o 2e a. te 
Goyy ~ (—=*) )d( 2 - (4) + 


5 Cac Ce Be 
Go,, ~ (A) )a(_S - —_ + 


eo, dexy TO, dev, ye Ce ; 


which simplifies to 


_ # glde,, (20, ,-(oyyt0,,)) + dey y(20,.-(o,,40,,)) 
+ de, ,(2o,, aie, NE + PNG afte en hte 
o74€,,) 
The term O45 dess may be expressed as 
o;7 des, = eo (oyyto,,) de, 4 


3 


2 | 2 
(ge yy" lo, #0, 7) dey +(Go,, P (9, to) de, 


3 3 


-\- 
oxy deyy FOV, de. + a5, des.) 
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Therefore 


“ Cee Parc 
ey) de. O.-. | 


Now the deviatoric stress increment do. § from equations (3-16), 


(3-17), and (3-19) is equal to 


de 
Peel dc. eeman ls k] 
do. G( ae Os. ~ ) 


= 26(ders - 855 Sep - 975 SRT eKI) (3-20) 
a > 


5 


which was obtained by recalling equations (3-6) which may be expressed 
in incremental index notation as 
Ge... =wGe seer O..: de 


ij ij ioe kk (3-21) 
3 





and substituting des 3 and dA into equation (3-7) and solving for do. 5 
equation (3-20) was derived. 


By definition the total stress increment do, « 1S 


IE 
do... = do.< + 6.. de _— 





do.% se G 6..de 


(3-22) 
Tj 3(1-2v) 1j kk 


Finally substituting (3-20) into (3-22) 


de 
2? a 6. ;¢ oy, 
do; 6 — E4470 ij <k “in ) (3-23) 


which may be expressed in matrix notation as 


{do} = ELD Jide} = 2(14v)6[D' Tide} (3-24) 





Mention should be made here of the particular notation for the elastic 
and plastic stress~strain equations. The elastic equation is expressed 
in the usual tensor form while the plastic equation is expressed in 


differential form. The plastic problem is point dependent, while the 
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elastic problem is not. The value of [o°} is a function of what point 
on the plastic stress-strain curve the calculation is being made while 
the value of ro is a constant anywhere on the elastic curve. 

The elastic and plastic equations are similar enough to allow the 


replacement of [oe 7 for a yielded material by cory. Explicitly ror 


























becomes 
]-v a oo - v Come (chy We! CuaG | GuaG 
ie x _ XX yy eee 22 OXY | XK Vg ene 
]-2y S | l-2v S ]-2v S S | S S 
| 
]- oe aoe ¢ : ¢ 
| YY “yh 2z)_ Fy xy |. “yy"yZ_ “yy!2x 
]-2v S 1-2 S | S S S 
2 ae. 
aes 
. T-2v S | 
ia ! | 
tv | 
ne 
| 
2 
Foes) Smee 
| 
Symmetry 











Le oF ~ Oyz°2x 
se ” i 2 (3-26) 
I+v ~ O_O ] o. 4 
ZZ oe 2x 
S : S 





The plastic stress-strain matrix [5 a unlike the elastic stress-strain 
matrix ro87, does not uncouple to give independent xz and yz behavior. 
The matrix roP is the non-linear, point dependent, plastic stress- 
strain matrix relating incremental stress and strain. 

The next section describes how roe and co" 7 fit into the scheme of 


the analysis to determine the failure point of the composite. 
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IV METHOD OF ANALYSIS 


“mes 


The objective of this analysis is to determine the failure point of 
a unidirectional fiber reinforced composite material subjected to longi- 
tudinal shear loading and to determine the effect that the matrix 


material has on the overall properties of the composite. 


A. THE NON-LINEAR PLASTIC PROBLEM 
The method for the plastic analysis described here was developed 
by Yamada, Yoshimura and Sakurai [16] for plane stress and modified for 
the longitudinal shear loading case. Briefly the analysis solves the 
incremental force-displacement equilibrium equation 
[KP] {dw} = {dF) (4-1) 
where cK] is the stiffness matrix, {dw} is the differential displacement 
in the Z-direction and {dF} is the differential force vector in the Z- 
direction. Equation (4-1) is obtained by finite element discretization 
of the displacement Ween fuk, this may be written : 
[KP (w)]{dw3 = {dF} (4-2) 
and solving for {dw} 
. {dw} = [KP (w)]7 {dF} (4-3) 
In the above problem the differential displacement is calculated on the 
basis of the preceding displacement, that is, the non-linear problem 
is solved by incremental theory, 
{aw,} = CK (we ,)7" {AF} (4-4) 
The stiffness matrix cry in the preceding discussion is a function 
of the displacement vecter {w} in the plastic analysis only. The stiff- 


ness matrix for the elastic analysis is a function of two independent 
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elastic constants and the equation {dw} = [K] {dF} may be integrated 

to give {w} = Key {F}. Hence the linear elastic problem can be done 
in one analysis while the non-linear plastic problem requires repeated 
incremental analysis. That is, the non-linear problem is approximated 


by a succession of linear problems. This is illustrated in figure 2. 





IN A (ELASTIC RANGE): F=kw 
IN B (PLASTIC RANGE): AF=kaw 


FIGURE 2. THE ELASTIC-PLASTIC PROBLEM 


— 
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The present discussion is limited to general terms, but a more 
detailed presentation will be related in Section 5. 

The problem considered in this analysis is the longitudinal shear 
loading case. The edge x=a 1S given a uniform unit displacement in the 
positive Z-direction as shown in figure 5. The resulting macrostresses 
Se and Sz and microstresses Oy5 and O72 are calculated by the finite 
element method using a linear strain triangle (6 degrees of freedom). 
Felippa [6] presents a derivation of the linear strain triangle for in- 


plane behavior and the derivation of the linear strain triangles for 


longitudinal shear behavior is presented in Appendix A. 


B. DESCRIPTION OF THE FILAMENT COMPOSITE 

The microscopic stress state in a composite material depends upon 
the arrangement of the filament within the matrix and the applied loading. 
To simplify the present analysis, we consider the special case in which 
the filaments are arranged in a doubly periodic rectangular array. This. 
array of filaments in the matrix may be divided into basic unit cells 
of dimensions 2ax2b, to facilitate ease of calculation. - This arrange- 
ment is illustrated in figure 3. The arrangement illustrated in figure 4 
with square cell a=b=0.0052 in., and filament radius c=0.004 in., results 
ina filament volume of approximately 50 percent of the total volume. 
The symmetry of the basic unit cell allows the analysis to consider only 
One quadrant of the basic unit cell as shown in figure 4. The quadrant 
considered has dimensions axb with a filament of radius ¢c. It is aiso 
assumed that the filament and matrix material are homogeneous and 1S0- 
tropic and that a perfect bond exists between the filaments and the 


matrix material, i.e., there is displacement continuity along the interface. 
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FIGURE 3. RECTANGULAR ARRAY OF FILAMENTS AND A BASIC UNIT. 





FIGURE 4. A QUADRANT OF THE BASIC UNIT. 
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The finite element model used is shown in figure 6. The model is 
divided into 92 elements and a sya of 215 nodes or degrees of freedom. 
The 92 element model was chosen over a 54 element model and a 200 element 
model. Previous analysis of Lin et al [11] indicated that the 54 element 
model did not yield results of sufficient accuracy, while the 200 element 
model gave very accurate results, but were not commensurate with the 


amount of computer time required. 


C. DERIVATION OF THE EQUILIBRIUM EQUATION 
In the elastic region the resulting strains obey the linear Hookian 
stress-strain relation, 


igh fl) fe 
6x1 6x6 6X1 Je 
where {co} is the column vector of stresses, cof is the elastic stress- 
strain matrix and ca is the column vector of elastic strains. The 
explicit form of the symmetric matrix ro87 is given by equation (3-1). 
The stiffness matrix [k] for the elastic elements may be derived from 


energy considerations. If 


dU = Oss d E. 


J J’ 
then the principle of virtual work shows that, 
[k] {w} = aU/. 6, (4-6) 


The particular expression for the stiffness matrix of the longitudinal 


Shear problem for the elastic problem is obtained as follows: 
Us [|[tortaerav = 2 | ltosTtesay (4-7) 
V V 


where we used the fact that stress is linearly proportional to strain, 


where the (om is the transpose of {o}. Using equation (4-5) and (4-7), 
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we obtain 


u a) | unseat | i (4-8) 
V 
The strain-displacement relation may be expressed as 
{e} = [B]iw} (4-9) 


where [B] is the coefficient matrix and {w} is the displacement vector 
in the Z-direction. Substituting this relation into the previous 
energy equation (4-8) | 

U = oo] [t2D 00°C t dA{u} (4-10) 


where t is the uniform thickness. From equation (4-10) the elastic 


stiffness matrix for an element becomes 


rk = {f te" to®7"rettan (4-11) 
A 


where A is the area of the triangular element. Similarly for the plastic 
element, the differential stress-strain relation is 


fdo} EPDM {de} (4-12) 
6x] 6x6 6x] 


where the explicit form of [oP is given by equation (3-25). In this 
case the element stiffness matrix is 

cePy =|] cei" co"y"centan (4-13) 

The stiffness matrix cKey for the elastic elements and ck for the 
plastic elements are assembled to form [K], the stiffness matrix of 
the whole body. The system stiffness matrix relates the incremental 
load {AF} to the incremental nodal displacements {Aw} for the elastic- 
plastic analysis 
{AF} = [K] {aw} | ala 

Equation (4-14) is basically how the plastic, non-linear problem is 


approached. The non-linear problem is linearized by taking very small 
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finite steps (increments) along the plastic portion of the stress-strain 
curve. In this way the equivalent stress for each element will not 
deviate during succeeding steps from the actual: stress-strain curve to 
any significant amount. The following discussion supplies a little 


more detail to this sketch. 


D. SOLUTION TECHNIQUE 

Using the incremental load technique the elastic-plastic problem is 
approached as follows. First the elastic problem is done in one step. 
For a test shear displacement (unity), solve [K]{w} = {F}. This equation 
holds for the elastic range only. Calculate the elastic displacements 


at the nodes and then the elastic strains, stresses, and equivalent 
E 


stresses 5° for each element. Determine the maximum o- and scale all 
elastic values using a scale factor of r = V/5 ay? where Y is the yield 
stress of the material resulting from a uniaxial tension test. This is 
done to bring the first element or elements to yield. Next calculate 

the fo’ and resulting rk for the yielded element or elements and re- 
place the previous elastic matrices roe and rk}, with the new matrices 
in the system stiffness matrix [K]. Next a test increment fin of the | 
shear displacement is chosen (where the superscript T inside the brackets 
dente Ge test increment). The strain increment be; 5! and the stress 
increments bo; are calculated at each element after solving for taw!}. 


Now scale all elements remaining in the elastic state by the minimum of 


r+ (r24a(ag T)° (y2-82))2 

a a ee as 4-15 
— (4-15) 

where 


ree (ag Omen esas! = (xe! )2 (4-16) 
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T P 
ACs; 


Ao = (2 Ao. (4-17) 

R is the scale factor that brings the next elastic element to yield 

(see Appendix B for derivation). Also, a is the present equivalent 
stress of the elastic element and ng! denotes the increment of co induced 
by faw' }. 

A fine point should be made clear before proceeding, ie and ae are 
not the same quantities. A significant difference does exist between 
the change in effective (or equivalent) stress fie!) and the effective 
change in stress (Ko'). The change in effective stress is the increment 
of equivalent stress (Gc) caused by an incremental increase in the loading. 
The effective change in stress is that portion of the stress increment 
that actually produces a change in the plasticity of an element. 

At this point in the procedure the nOdaUmen ep lacencnren strain and 
stress increments are multiplied by the minimum value of R and addea to 
the existing values. Next the equivalent stress o, given by equations 
(2-9) and (3-11) and the incremental equivalent plastic strain, ee are 
calculated for each element. 

= O45 Aeiy 
ait at Linea 


If the value of the ne’ is positive, the analysis returns to the be- 


(4-18) 


ginning of the next plastic cycle by calculating co" 7 and ck] for the 
next yielded element. If wa os negative the computation is stopped. 

This is included here because large plastic strain are possible without 
any practical increase in the load after the plastic region has expanded 
tremendously. This gives an indication that the analysis is very close to 


the collapse load of limit analysis. The problem may terminate in a very 
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natural mode when all the cycles are complete or when the incremental} 
equivalent plastic strain, hes goes negative. These two modes of failure 
are illustrated in figure 7. The next section describes the actual 


procedures in much more detail. 


<a 


FAILURE 
POINT 





FIGURE 7. MODES OF FAILURE 
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V PROCEDURES OF CALCULATION 


This section provides a much more detailed explanation of the pro- 
cedures of calculation than does the Method of Analysis section. The 
procedures of calculation of this analysis was suggested by Yamada et al 
[16] and modified for the longitudinal shear loading problem. The analysis 
was constructed to be able to perform calculations for both the xz and 
yz shear loading problems, but will only be described for the xz shear 
problem. The xz shear problem was the one actually considered by this 
investigation. The procedure follows. 

For the elastic analysis (Kw=F), 

1. Apply a uniform unit displacement to the x=a face, and calculate 
the resulting displacements, strain, and stresses for the elastic be- 
havior (system stiffness uses elastic stress-strain relations here). 
Figure 5 illustrates the uniform unit displacement given to the basic 
unit cell (first quadrant only). 

2. Scale the elastic values so that the element with the maximum 
equivalent stress (') is at initial yield (Y), i.e., adjustment of 
input displacement so that a single point is at yield. The scale factor 
re 1S | 


es 
oe Y/=7 (5-1) 
max 


and the nodal displacement vector (we } and the nodal load vector (ay at 
the initial yield point are given by 
we} = rFqw) | (5-2) 


Fey = rege!) mn 5-3) 
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This terminates the elastic analysis. For the plastic analysis that 
follows we proceed in cycles or increments. The procedure for the plastic 
analysis is one cycle that is repeated over and over until one point or 
element of the composite has failed (termination of the problem may end 
in other ways, see section four and actual computer program). The loading 
for the plastic problem may mean either displacement or force loading, 
but for this particular study the displacement loading was employed. 
Displacement loading is used because we are interested in achieving the 
shear fetoemn tenis tate. Thus we use equation (4-4) for this problem 
({dw.} = [K(w,_,)] 7 (dF). 

3. Calculate De and rk for the post yield elements; usually 
only one element will reach yield in the first cycle. However since 
elements close to yield will be carried into yield during the next incre- 
ment of displacement, any element that is within .995 of yield is con- 
sidered to be at yield. The value .995 is completely arbitrary; it 
depends on the accuracy that the user wishes to obtain. 

4, Modify the system stiffness matrix to reflect the stiffness of 
the elements that have yielded. 

5. Choose a test increment aw! of displacement. 

aa Solve tar! } = [K(w. _,)J dw," } and then calculate the strain 
increment {Ac} and the stress increment {Ac}. 

7. Calculate the scale factor R (equation 4-15) for each remaining 
elastic element to reach yield. 

8. Determine the minimum R. The displacement increment, R «(OW 
is sufficient to yield that element having the minimum R and place it in 


the plastic range. 
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9. Multiply the strain increment and stress increment calculated 
in step six by Rede and add the resulting del (A) quantities at the 
nodes to the previous nodal values, i.e., {wh s={w},  +{Aw}. Store the 
results. 

10. Calculate and store the new equivalent stress for each element. 

11. Calculate the increment of plastic strain for each post yield 
element. The stiffness of each plastic element constantly changes as the 


element progresses into the plastic range. 


-— ee AE T 
Pe (4-18) 
o(1+H"/ 4.) 
For this analysis —— (o._ Ae._to_. Ac_.) 
1? VZ = YZ 42% 2 ZX (5-4) 
o(1tH"/,0) 


12. Check that ae positive. If nes 4c positive return to step 
three. If fee Ge negative stop computation. This is included because 
large plastic strains are possible without any practical increase in the 
load after the plastic region has expanded tremendously. This is one 
way of ending the problem, the other way of course is a natural ending. 
A natural ending occurs when the number of plastic cycles requested are 
run or where some of the other checks built into the program terminate 


the problem (see computer program in Appendix C). 
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VI RESULTS, CONCLUSIONS AND RECOMMENDATIONS 


A. RESULTS 

The computer program for this elastic-plastic analysis of the 
longitudinal reinforced composite loaded in the xz shear mode was run 
until some point in the matrix reached its ultimate strength. The 
properties for the matrix material and for the boron filament used in 
the analysis are given in table 1. 

In the analysis the first point (element) to yield was calculated 
before the ultimate breaking point was located. The initial yield area 
for boron-epoxy and boron~aluminum composites occurred at element 50 
as enon in figures 8 and 9 respectively. Bordering the filament~matrix 
interface, this point (element 50) Ties in the matrix material on the 
lower horizontal face of the model. The filament point in both com- 
posites that obtains the highest equivalent stress at the time of initial 
yield is element 32, which is the filament element adjacent to matrix 


yield element 50, (see figure 6). 


Table 1. Material Properties 


Property. Boron Aluminum Epoxy* 
Young's Modulus 60x10° psi 10x10° psi .5x10° psi 
Poisson's Ratio 0.20 0.30 0.31 
Elastic limit 480,000 psi 34,600 psi 3000 psi 
Ultimate Strength 480,000 psi 42,500 psi 4200 psi 

N = 38.3 Or Wi 
B = 46,300 psi 5088 psi 


* Tensile Properties for Epoxy 


4) 





When the analysis finally reaches the breaking point for both com- 
posites, some interesting results have occurred. In the boron-epoxy 
composite element 50 becomes the first matrix element to fail and 
element 32 is still the highest stressed filament element, but in the 
boron-aluminum composite the initial yield element 50 is overtaken by 
element 70 to become the first failure element and element 32 is over- 
taken by element 31 to become the highest stressed filament element at 
the time of failure. The progression from initial yield to failure for 
the boron-epoxy and boron-aluminum composites is shown in figures 8 and 9 
respectively. 

The results of the initial yield analysis are presented in table 2. 
These are the results associated with the elastic limit analysis. The 
boron-aluminum composite has a macro shear stress Sy 77135372 Sima tees 
elastic limit; and the boron-epoxy composite has a macro stress of 


oa O48 psi at its elastic limit. 


Table 2. Results at Initial Yield 


Boron-aluminum Boron-epoxy 
hee, 19,953 psi 1729 psi 
Snax so 7 34,600 psi 3000 psi 
Snaxa, 35,288 psi 3070 psi 
Se 13.372 en 1044 psi 
Svz 240 psi 22 psi 


We note that the maximum stresses in the filament are but a smal]? fraction 
of the elastic limit stress (480,000 psi) of the filament material, that 


is, the filament is not being used efficiently. Thus, designing to the 


42 








FAILURE 


Figure 8. PROPAGATION OF PLASTIC ENCLAVE FOR BORON-EPOXY 
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6. FAILURE 


9. PROPAGATION OF PLASTIC ENCLAVE FOR BORON-ALUMINUM 
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elastic limit of each, the boron-aluminum composite can sustain a shear 
load 1300 times greater than the boron-epoxy composite. 

Table 3 contains the post yield results. Continuing the loading 
into the plastic region, the more ductile aluminum matrix allows the 
boron~aluminum composites’ longitudinal shear loading, 25 to increase 
80 percent (from 13,372 psi to 24,140 psi), while the boron-epoxy 
composites’ ae increases 50 percent above its elastic limit (from 1044 psi 
to 1539 psi). These increases in longitudinal shear load correspond to 
a 25 percent increase in the equivalent stress of boron-aluminum (the 
equivalent stress in the aluminum element 50 goes from 34,600 psi to 
42,500 psi), while boron-epoxy's equivalent stress increases 40 percent 
(s increases to 4200 psi from 3000 psi). Clearly the boron-aluminum 
offers the designer a better shear load carrying ability than does the 


boron-epoxy composite when the design is allowed to reach failure. 


Table 3. Results at Ultimate Strength 


Boron-aluminum Boron-epoxy 
x76 24,444 psi 2427 psi 
Sy 70, 24,563 psi - 

OnaX «9 42,372 psi 4212 psi 
Cmax 42,547 psi ~ 
OmaX a0 51,097 psi 4356 psi 
Oak D2, /S emp SH ~ 
ae 24,140 psi 1539 psi 
Syz 692 psi . 33 psi 
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An all aluminum unit cell was subjected to the longitudinal shear 
loading to verify the analysis’ reliability. A Ke macro Shear of 
20,000 psi was obtained for an equivalent micro stress, o of 34,600 psi 
at initial yield. The homogeneous cell had a uniform microstress state 
throughout, with o,, = s//3 = 20,000 psi and oy, = 0. The ultimate 
macro shear stress ae is 24,500 psi with a uniform equivalent micro 
Stress o = 42,500 psi, resulting from as 24,000 psi and G7 = 0 psi. 
We note in the case of a homogeneous cell subject to pure xz Shear, 
there results only xz behavior, that is, in this special case post yield 
behavior does not give coupling between xz and yz behavior. This is 
evidenced by Ov = Q psi. 

The composite results were further verified by constructing an 
equivalent stress-strain curve (universal stress-strain curve) and then 
by plotting the values of equivalent stress and strain from the computer 
program on the same figure to show how close the analysis came to the 
input stress-strain curve. The elastic portion of the curve was con- 
Structed using the linear elastic relation, o = Ee, while the plastic 
P (SN 


portion of the curve was constructed using the power law, € = B 


The portion of the boron-aluminum equivalent stress-strain curve where 


the slope becomes constant (i.e., for & > 40,500 psi) was plotted using 


the relation, H* = The power law values, N and B were calculated 


do 
using the data from’&n original uniaxial stress-strain curve for the 
material. The determination of N and B for the aluminum was easily ob- 
tained, but N and B for the epoxy curve became extremely difficult to 
calculate due to the almost linear stress-strain curve for epoxy in the 


plastic region. Figures 10 and 11 show the equivalent stress versus 


total strain curves for the boron-epoxy and the boron-aluminum 


~ 
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respectively, while figures 12 and 13 show the equivalent stress versus 
the equivalent plastic strain. The close proximity of the plotted 
values to the curves verify small enough incremental steps were taken 


during the non-linear plastic analysis. 


B. CONCLUSIONS 

In this analysis the objective has been to develop and use a program 
which could determine the failure point of a unidirectional reinforced 
composite material subjected to longitudinal shear loading using the 
finite element method and to make a comparison of two composites. 

The failure point was determined for the two composites and occurred 
at element 50 for the boron-epoxy composite and at element 70 for the 
boron-aluminum composite. A greater portion of the aluminum matrix 
yielded before failure than for the epoxy matrix. ThisS occurrence came 
as a result of the more ductile aluminum matrix. 

The appearance of the longitudinal macro shear Sy is smal] compared 
to the longitudinal shear Sy? but must be accounted for because the 
Ove microstresses were comparable to the one microstresses for the inter- 
face elements (values of 50 percent or greater were reached). Since the 
significant Ove stresses occur in the interface elements when the model 
was Tosced in — Shear ,xz- yz coupling may be explained as a result of 
the circular filament cross section. If the filament cross section had 
been rectangular the C7 terms and hence Syz term would have been less 
likely to appear. 

Considering all the results the boron-aluminum composite definitely 
Seems to possess a better longitudinal shear load carrying ability than 


does the boron-epoxy composite. But the all aluminum cell gives the 


Same results. Hence it would not be economically wise to use a boron- 
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aluminum laminar composite in place of a homogeneous aluminum plate if 
in-plane Shear loading was the only load applied. The boron-aluminum 
composite may not fare to well when compared to an all aluminum material 
in the longitudinal shear mode, but from the previous analysis of Lin 

et al [10] we see that the all aluminum material would not approach the 
boron-aluminum composite when loaded in other modes (i.e., transverse and 


longitudinal loadings). 


C. RECOMMENDATIONS 

The logical extension of this analysis would be to combine it with 
the other loading problems studied by Lin et al [10] to construct a fail- 
ure surface. Also the investigation of different shaped cross sections 
for the filament should be considered. These suggestions combined with 


a parametric analysis could lead to other interesting thesis topics. 
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APPENDIX A 


CALCULATION OF THE STIFFNESS MATRIX 
USING FINITE ELEMENTS 


Section three derived the stress-strain matrix for the plastic analysis 
and stated the stress-strain matrix for the elastic analysis. This sec- 
tion will show how these stress-strain matrices are expressed in finite 
element form so that they may be written in computer language. 

In this analysis, the longitudinal shear loading problem is associated 
with the macroscopic stresses Sy and Suz and the displacement u=v=0, 

=w(x,y). The finite element used in this analysis is the linear strain 
triangle (LST). For the longitudinal shear element the elastic strain 


energy 1S given by 


U = AIT, se, 5a (A-1) 
y 


which becomes for this problem 


U = 5 |i] (20 e125 se0_ dV 
a. | MZ X Ze Y ZY Z 


| V 
- 4 
“on all (2% x2*yz yz) (A-2) 
V | 
where oe and ov are the microstresses and spre and Yyz are the shear 


strains associated with the problem. Assuming a unit thickness 


t | 
U=£}/6,,0 "ys, V xz dA= 
A yz . 
t {Je xz{!J¥xz{dact t HF tox" ty dan (A-3 
A Dire y 
Al-yz yz A 
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The minimization of the strain energy U with respect to the displace- 
ment vecter {w} gives the stiffness matrix of an element. The derivation 
of the element stiffness matrix follows. 


In the elastic elements, 


o. = Gy. and o Gy 


XZ XZ yz yz 
6} GO] ly 
y2 ¥2 
Substituting into equation (A-3) 
aay [oe ]iy}dA (A-5) 


Before the derivation may c opr the array {y} must be defined. 


| Y tC, GS b 0 0 O i 
re ay ee es XZy 


: 
yz 0 0 0 by be 83 "x2, 


’ 
= C XZ 
2x 6 "yz, 7 


where Y, are triangular coordinates. 
The ees and V2 ° go from one to three for the corner nodal points 
of a linear strain triangle. Now the nodal shear strains may be ex- 


pressed in terms of the nodal displacements 
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: ww _ 9 ‘ 
fy, } = = — CN NONZNENENGS [| 4- 
“iat ox ox < > j We 
Wy 
=<N NON NNO NN A-7 
<i Xo Xg Xy Xs xe> j We (A-7) 
and vy 
tyz.) = cA = & A NaNstaNlsNg> | 
L ay ay 1 [We 
Wi 
=<N NON N NN ; (A-8) 
Dy Nyollyaly ly. Y6> , 7 
where the Ny and NY are the first partial of the shape functions with 


7 1 
respect to x and y respectively. Evaluating equations (A-7) and (A-8) 


at each of the six nodes, 


"XZ1 06 . Qa Sh, a an ed 


e*ee 


*@¢¢ 


= Gi N ' .-10 
Tyzq 6 ee a (A ) 


Combining results and putting equations (A-9) and (A-10) in matrix form 


{yo} = L2x] {w} 


A-1] 
6X1 a eG 
fy 7} = [zy] {w} (A-12) 
. 6x1 6x6 6x1 
Substituting equations (A-11) and A-12) into (A-6) we have 
ex] 
{y} = [é] {w} (A-13) 
[zy] 
= Ce] £2] ow) | 
Now using the equation (A-13) and equation (A-5) we finally have 
u -£) bo" t2I ter WN Leatedimen (A-14) 


Letting [B] = [e]fz] we have 
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u = rrew3' £8)" (oFy ced widA (A-15) 
A 


Minimizing U with respect to {w} and integrating 
[k“] = [B]'[D"ILBItA (A-16) 
apace is identical to equation (4-11) 
For the plastic elements the stiffness matrix becomes more involved. 
The stress-strain relatives are a function of the loading (see equation 
(3-26)). The only difference is that ror] takes the place of roe in the 
plastic stiffness matrix, thus 


cee] = cel! co ypeyta. (A-17) 
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APPENDIX B 


SOME DERIVATIONS OF EQUATIONS 


In this appendix will be presented some important derivations necessary 
to the understanding of the elastic-plastic analysis. There are two 
proofs of the scale factor R used in section four and five. The first 
proof is geometrical in nature, while the second is an algebraic deriv- 
ation. The derivation of the gamma, which is necessary for the scale 
factor R derivation is also included. Finally it is shown that the con- 


stant C of equation (2-18) equals 3/2. 


A. GEOMETRIC DERIVATION OF EQUATION (4-15)- — 
Here we obtain equation (4-15) through a geometric viewpoint. The 


basis of the discussion 1s contained in figure 14. In figure 14 o is the 


present equivalent stress of the elastic element and ig denotes the in- 


crement of co induced by the displacement increment uae. o and nee are 


od 


represented respectively by OP and PS, while Aco Ve PR. For the longi- 


, eee a 
tudinal shear problem ae (see equation (B-7)) and (Ao Tye (see 


equation (p-13)) determines the load increment sufficient to just cause 
yield in each elastic element. 


Wane 2 Ae 


al | ol 


Using the cosine law, 


Let x = o2+(Ac 
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FIGURE 14. DETERMINATION OF THE LOAD INCREMENT TO 
JUST YIELD THE ELASTIC ELEMENT. 


and 
COSn =x = ARCCOSx SINa = Y 1-x2 


Using the sine law, 


lS 
© 


pos Us 6 = 180 - (a+60) 


SINa  SIN®  SINO 


0O=Y OP=3 


Ses ae ee 
J/1-x2 SING 1-x2 — SINO 
= Se? 
sino = © (VIX) Lg 
ti 
© = ARCSINZ 


08 


IS) (ae: 
al [sl 








¢ = 180 - (ARCCOSX + ARCSINZ) 


See QC‘ 
/1-x2 SIN(180- (ARCCOSX + ARCSINZ)) SIN(ARCCOSX a ARCSINZ) 


LET g = ARCCOSX and Y = ARCSINZ 


Ene 
SINBCOSy + COSgSINy 


PQ ale 
( 


iH 


P 
[ 


22, 
7) 
I 


Y{/1-z2 + ar (B-1) 
-X 


v4 1-(8)(1-x) 4 x (SN(1-x° yy 
J 1-x2 
2 2 2 


Y - 6 (1-x ) + xa 





iN 


il 
PO = {[(v2-82)4(a3 1)?+x27/4(03 1)7)? 4 


2 9 2 _ 
(3° )+(as_ ') -(stas!) y/a(as *) (B-2) 
Expanding terms and letting 
_ 2 2 
T = (Ao yy TENS - ie (B-3) 
PQ becomes 
eT + 1} 
PO = [{(Y2-62)4(Ao ')2+r2}74+7]/ 206 (B-4) 
PQ T ae 
R= ——= [{(Y2-c2)4(f5  )4+r2}4+r) J 
PR a 


aume ye (B-5) 


B. ALGEBRAIC DERIVATION OF EQUATION (4-15) 


The next derivation of R may be expressed in general terms but jis 


instead derived with the particular problem of longitudinal shear loading 


a mince 


oo 








When g+taG = Y yield occurs, Aoy is produced by the incremental test 
displacement awl. 


ee Z Z - 
O 3Lo, toy ] (B-6) 


(s+a5')2 = 3[ (6 thoy, )?+(oy,*ho Ty2q (B-7) 


x yz 


6*+26 ages is elo Eco ACs 


+ (Ao ere 2425 _ Ao i One yy 
XZ VZ dae VZ aaVZ yee 


Also when Gaye = Y2 element is at yield: 


‘F =e | 
Za = 2 ae 2 2 
Y 3[o,, +2Ro, ING +R (A5,,, ) 


ae | a il 
Z Z 2 a 
+ Ve BCAoh As, +R (A5\,5 )2] (B-8) 


hee, 
Y O +3[2R(o, AG, to, Evy.) 


1 INEUIG +ao, J") 


7 ee |, - T 
T = -6[o, AG, tO, AG ] (B-9) 
—T a Seah 
(Ao )4 = 3L (as, )*+(d5., 2] (B-10) 
“Making the necessary substitutions, 
(Ko" )2R2-rR+52-Y2 = 0 
R = r+y 7244 (ho! )2(y2-32) (B-11) 
ATE 
C. DERIVATION OF EQUATION (4-16) FORT 
From the original postulate 
ve (re epee insta (B-12) 
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——T 


AM 
Ago = (Sho aD 


2. 0. 2 


ij 1j 
which equals PR in figure 14. The term ice becomes upon expansion 


— | a, -T2 “72 “72 st 2 sali 2 
Ao = {5 (do,,, vai vcr *2( (Ao, y ) eV ) 


i 
2 


+ (do, ')2))3 (B-13) 


The equivalent stress vector is 


a 1) +2 22 2( 22 ; 
4<G2 = 6g + +o +2\(o to to 
3 XX oyy ZZ Xy yz XZ 


(Staal yeb= 62426 Ao! +(ag!)2 (B-14) 


putting co + fel in the form of the equivalent stress vector 


pea ee ey 2 oe "eo eee 
3 (G+AG ) (,  +AG, ) TA OE ae ) + 


7 = Pee 7 me ee 
Gee thc. ) #206, tA5, y } 


= OE a 
: ‘Sy A055 ) J 


7aas +2(03")? : 3 (3+as' )?-26? (B-15) 


Now the outcome is subtracted from (eye to give r. For this analysis 
T becomes 
Te itenucline 2 


Wee eee ee | 
= O15, y Ao, to. Aco _ al 


y “yz “yz Eel) 


D. DERIVATION OF C IN EQUATION (2-18) 
The proof that constant "C" of equation (2-18) is simple and straight- 
forward. Equation (2-19) for a tension test and principal strains 


becomes 
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ede de. jee 
faq =f y (de, de.) 
p aa ep 
[254 (dey? )24(-M2_y24(—861_y2yz = J Oe ae 
Z 2 


and from equation (3-11) the equivalent stress becomes 
o = [3(o;; 65; 12 = [Boy “tog toy 12 (B-18) 
where ard and 5S are the principal deviatoric stresses. Equation 
(B-18) simplifies to 
Onis J ((804)2+ (hoy )2#(-40y )232 = o1 
and for tension i = £3. Now substituting into equation (2-17) 


a 


G 


a oe 


o,? 


P 
: 


ie) 
2 


G 


wir 
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APPENDIX C 


THE COMPUTER PROGRAM 


The computer program performs an elastic-plastic analysis for a uni- 
directional reinforced composite material subjected to longitudinal shear 
loading. The main sections of the program are listed and described below. 

Main - Constructs the elastic stiffness matrix and modifies it according 
to loading and boundary conditions. It calls the equation solver and then 
scales the elements in the elastic and plastic analysis and controls 
further calling sequences. 

Post ~ Constructs the plastic stiffness matrices for the elements in 
the plastic range. 

Caltau - Calculates element strains and stresses: c=Bw, o=DEe and 
o=D' ec. 

Bandec - Solves banded matrix problems (Kw=F, or K(Aw)=AF). 

MULI - Multiplies a symmetric banded matrix with a vector. 

| This portion of the discussion of the computer program is devoted to 
assumptions made in writing the program and are 1isted here for 
convenience. — 

Ve The properties of the matrix material may be different in com- 
pression and in tension, YMC and YMT respectively. This assumption 
explains the dilatation variable SII‘'s inclusion in the program. SII 
was the variable that determined if a matrix element was in compression 
or tension. It was included to preserve the generality of the program 


but does not enter the analysis in the particular case of longitudinal 


Shear loading. 
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2. If an element is within a particular fraction (called TEST) of 
the yield value, it is included as a yielded element. In this analysis, 
TEST was chosen as 0.995. Therefore for calculations any eljement which 
is at .995 times the yield value is at yield (TEST = .995). After an 
increment it can be determined whether an element which was close to 
yield and assumed to yield, actually did. In other words, the size of 
the incremental step and TEST should form a consistent pair. 

3. For purposes of scaling the load to take the critical element to 
yield, the scale factor Beas is restricted to be less than some maximum 
(RYMIN). This is imposed to limit the size of an incremental step taken 
along the stress-strain curve. If too large size steps are taken, the 
elements in the plastic region will provide inaccurate results (see 
figure 2). Also at a certain point in the stress-strain curve larger ane 
values will be possible which would reduce the calculation time. These 
maximum values for Keay along the stress-strain curve must be altered 
according to the materials being used. They are also a result of experi- 
ence in working with a certain material to produce better results. 

4, Also in scaling the load to take the critical element to yield 
the filament elements will be ignored. The yield of a filament is much 
higher than the yield value for a matrix element. 

5. For purposes of computer efficiency the elements of the stiffness 
matrix were derived to facilitate their use in the construction of the 
stiffness matrix in the main program and the strain column matrix in the 
subroutine Caltau. Therefore the zx and zy variable arrays listed are 
the first partial of the shape functions with respect to x and y 


respectively. 
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Cah at 1 ei Sat 1 
aX dy 


Pei > 2 ; ai ip 5s 2 
” ae [zy] = x 
x | 
ali Se 3 wm ie 3 
aX AG 1 See 


Cee FOC a dc/aee a two point curve fitting scheme is used to 
approximate the plastic portion of the stress-strain curve. For example 


assuming the plastic strain-stress relation can be expressed as the power 


N-1 
= B 
B= 8 3) 
P 
€ 


The values of B and N may vary for different material and for tension 


ee ayn 
low z = (=) 


= 


and compression of the matrix. The calculation of N and B results from 
solving the equation a (3) for two points on the plastic portion of the 
curve. The selection of the two points is important and should approxi- 
mately divide the non-linear portion equally. The points are chosen in 
this manner so that large errors do not result. This method is an accur- 
ate and simple procedure. 

Other curve fitting expressions giving more accurate results may be 
used. It is felt that more accurate results are not necessary and do not 
yield that significant of an improvement in the results for the amount of 
time that would have to be expended. 


7. Data input is another important aspect of programming. The only 


possible trouble area exists in numbering the nodes of a triangular 
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element. The usual convention is followed by numbering the corner nodes 
first and then the intermediate nodes in a counter-clockwise fashion. 


For example, the numbering sequence should appear as: 
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